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ABSTRACT 

We report on the discovery of a 590 days long term periodicity in the hard X-ray component 
of the microquasar GRS 1915+105 found from a comprehensive study of more than four years of 
RXTE observations. The periodicity is also observed in the hard X-ray flux observed by BATSE, 
and in the radio flux as seen with the Green Bank Interferometer and the Ryle Telescope. We 
discuss various possible explanations, including the precession of a radiation induced warped 
accretion disk. 

Subject headings: Physical data and processes: accretion, accretion disk — X-rays: binaries: — X-rays: 
individual (GRS 1915+105) 



1. Introduction 

The galactic X-ray binary system GRS 1915+105 
is the most prominent microquasar showing dra- 
matic variability in its light curve (Greiner et al. 
1996), quasi-periodic oscillations, phase lags and 
coherence behavior (Morgan et al. 1997, Muno 
et al. 2001). Microquasars are thought to be down- 
scaled analogs to quasars harboring a stellar mass 
black hole which accretes matter from a com- 
panion star. They show superluminal ejections 
(Mirabel & Rodriguez 1994) making the sources 
potential laboratories for studying accretion and 
relativistic jets near black holes (Mirabel et al. 
1992). GRS 1915+105 is the most energetic object 
known in our galaxy with a luminosity of ~ 5 • 10 39 
erg/s in the high state and ~ 10 38 erg/s in the low 
state. The binary system, located at a distance of 
~12 kpc, contains a 1M late type giant (Greiner 
etal. 2001a) in a 33.5 days orbit around a 14 M 
black hole (Greiner etal. 2001b), making it the 
most massive stellar black hole known. 
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Originally discovered by GRANAT (Castro- 
Tirado etal. 1992) GRS 1915+105 has been ex- 
tensively monitored by the RXTE satellite since 
1996, and a number of investigations of these data 
have been published (e.g. Greiner et al.1996, Mor- 
gan etal. 1997, Belloni etal. 1997 & 2000, Muno 
etal. 1999 & 2001, Rau & Greiner 2003; hereafter 
RG03). 

In a recent paper (RG03) we presented a com- 
prehensive study of the X-ray spectral behavior 
of GRS 1915+105 in the x-state (Belloni etal. 
2000) for over nearly four years of observation 
with RXTE. Within this study, a remarkable long 
term periodicity was found which was mentioned 
briefly by Kuulkers etal. (1997). They found a 
19 month long term variability comparing the X- 
ray outbursts observed with BASTE and RXTE 
data from August 1992, March 1994, October 1995 
and May 1997. This is consistent with our pre- 
sented result. 

A number of X-ray binaries, including LMXB 
as well as HMXB, exhibit long term variability 
in their X-ray flux with periods remarkably longer 
then their orbital period. Tananbaum et al. (1972) 
found a 35 days period in Her X-l which they ex- 
plained by the precession of a warped accretion 
disk in the system. The same process is used to 
describe the long term variations in SMC X-l (Wo- 
jdowski et al. 1998), LMC X-3 (Cowley et al. 1991) 



and LMC X-4 (Heemskerk & van Paradijs 1989). 

In this letter we discuss the detection of a new 
long term periodicity in GRS 1915+105. 

2. Data reduction and analysis 

The analysis presented here is based on the data 
reduction described in full detail in RG03. We in- 
vestigated public RXTE data from GRS 1915+105 
from November 1996 to September 2000 obtained 
from the HE AS ARC. We analyzed data of the 
PCU0 detector of the Proportional Counter Ar- 
ray (PCA) from 3-25 keV and data of the High 
Energy X-ray Timing Experiment (HEXTE) clus- 
ter from 20-190 keV for 139 %-state observations 
from 89 different days. x _ states correspond to the 
low/hard state of GRS 1915+105, exhibiting rela- 
tively low flux from the accretion disk and showing 
a hard power law tail. These states are character- 
ized by the lack of obvious variations in light curve 
and spectrum and are connected by continuous ra- 
dio emission of varying strength. 

For the spectral fitting in XSPEC 11.0 (Arnaud 
1996) we used a model consisting of (i) photo- 
electric absorption (WABS), (ii) a spectrum from 
an accretion disk consisting of multiple blackbody 
components (DISKBB) and (iii) a power law spec- 
trum reflected from an ionized relativistic accre- 
tion disk (REFSCH; Fabian et al. 1989, Magdziarz 
k Zdziarski 1995). 

The analysis of the X-ray model parameters 
for all x-statc observations of GRS 1915+105 re- 
veals a remarkably periodic behavior of the power 
law slope, T, which dominates the hard X-ray 
spectrum. It shows a nearly sinusoidal variation 
(Fig. 1, see also RG03) on a time scale much longer 
then the orbital period, P or b- The recently found 
correlation of the power law slope and the radio 
flux, Fr, in the x-states (RG03) suggests that one 
should also search for a long term periodicity in 
the radio data of the Green Bank Interferometer 
(GBI; 2.25 and 8.3 GHz) and the Ryle Telescope 
(RT; 15 GHz). 

For the determination of a cycle duration an 
analysis of variance for several model parame- 
ters and X-ray and radio fluxes was performed 
using the Fisher-Snedecor distribution function 
(Schwarzenberg-Czerny 1989). The results for T, 
the 20-200 keV X-ray CGRO/BATSE flux and the 
8.3 GHz RT flux are given in the periodogram in 
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Fig. 1.— 1.5-12 keV RXTE All Sky Monitor 
(ASM) count rate (top panel), 15 GHz RT (2nd 
from top) and 2.25 GHz GBI flux (3rd) from 
November 1996 to September 2000. (GBI was 
shutdown some time before JD ^2^50^00 and fol- 
lowing JD ^24-51660.) The lower panels show the 
power law slope, Y, (4rd) and the accretion disk 
temperature, Tf,b, (5th) for all observed x-states of 
GRS 1915+105. Error bars are la for each pa- 
rameter of interest. 
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Fig. 2. It shows a maximum at a period, Pi ong , 
of 590±40(FWHM)days. T and F R show sev- 
eral identical local maxima around 200, 300 and 
400 days which are caused by the three major ra- 
dio outbursts (see right panel of Fig. 3) during 
a 590 days period. This local maxima result from 
the time separation of the outbursts at phase —0.6 
and -0.85 (200 days), -0.2 and -0.6 (300 days) 
and —0.2 and —0.85 (400 days). These local max- 
ima in the analysis of the variance of Gamma 
arise from a combination of two factors. First the 
chi-states are usually connected with radio emis- 
sion and second the chi-states have non-stochastic 
distribution within the RXTE observations. The 
later factor results in an uneven coverage of data 
points within the time period. 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 — : 




EDO 400 600 BOO 1000 

Period [days] 



Fig. 2. — Analysis of variance of the 8.3 GHz GBI 
radio flux (top curve), the 20-200 keV BATSE X- 
ray flux (middle) and the power law slope (bottom) 
for JD 2450300-245 1800. The y-axis is in units of 
the test statistic of the Fisher- Snedecor distribu- 
tion function. The test statistic has larger values 
for the radio flux due to the higher number of GBI 
observations compared to BATSE and RXTE. A 
maximum at 590±40 days is clearly seen in each 
of the three graphs. 

Fig. 3 shows the power law slope and the GBI 
radio flux at 2.25 GHz displayed in segments of 
590 days each. These light curves demonstrate, 
that 590 days is the preferred period rather than 
200,300 or 400 days. Even though the coverage 
with GBI data is shorter than with RXTE ob- 
servations the periodicity is clearly visible in the 
radio outbursts. Both the broad outburst with ex- 
ponential cutoff at phase —0.2 and the shorter and 
steeper outbursts repeat periodically. 
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Fig. 3. — Power law slope, T, of all analyzed 
X~state observations (left) and of the 2.25 GHz 
radio flux from GBI (right) from JD 2450350- 
2450940, JD 2450940-2451530 and JD 2451530- 
2452120 (top to bottom panels) . Each panel covers 
a 590 days period and JD=2450225 was chosen as 
the zero point for plotting reasons only. The solid 
lines visualize the best fitting sinusoidal function 
(left) and illustrate the repeating outbursts in the 
radio flux (right). 

The soft X-ray component did not show any pe- 
riodicity in the %-states. This is true for the ASM 
count rate and hardness ratio as well as for the 
fit parameters of the accretion disk, i.e. effective 
temperature and the inner disk radius. 
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Fig. 4.— BATSE 20-200 keV flux from JD 
2448580-2452120 with a time binning of 57 days. 
The dotted lines devide the light curve into 590 day 
segments. 
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The 590 days period in the hard X-rays is also 
seen in the 20-200 keV BATSE earth occultation 
fluxes (Tab.l). However, the long term period- 
icity seems to be variable on larger time scale. 
Fig. 4 shows the 20-200 keV light curve with a 
52 day (orbital precession period of CGRO) bin- 
ning. This binning ensures the removal of any pe- 
riodicity arising from the precession of the satel- 
lite orbit. Although the 590 day period is also 
visible when analyzing the data from JD 2448300 
to JD 2450300 (before RXTE observations where 
performed, Tab. 1) the structures in the light 
curve show remarkable differences and a somewhat 
shorter period of 545±25 days is found when ana- 
lyzing all available BATSE data from JD 2448300 
to JD 2451800 together. A first hint for this vari- 
ation is seen when fitting the power law slope be- 
havior with a simple sinusoidal function with a 
590 days period (solid line in left panel of Fig. 3). 
The minimum in he power law slope of the long du- 
ration x-state from JD 2450400-2450600 lies out- 
side the minimum of the sinusoidal function (Fig. 3 
top). Also, an additional maximum is seen around 
phase=0.6, particularly in the bottom panel. The 
inpersistent behavior is not surprising, and prob- 
ably reflects the known transient behavior of GRS 
1915+105. 

Table 1: Results of the analysis of variance for 
several model parameters and energy ranges for 
the analysis over the time span of JD 2450300- 
245 1800 with exceptions explicitly stated other- 
wise, (n.d. = periodicity not detected) 



data 


days 


Tbb 


n.d. 


r 


590±40 


R 


590±40 


2.25 GHz GBI flux 


600±40 


8.3 GHz GBI flux 


590±40 


15 GHz RT flux 


600±40 


GBI spectral index 


590±40 


1.5-12 keV ASM cts/s 


n.d. 


20-200 keV BATSE flux 


600±40 


20-200 keV BATSE flux (jd 2448300-2450300) 


590±40 


20-200 keV BATSE flux (jd 2448300-2451800) 


545±25 



3. Interpretation 

The observation of the 590 days periodicity in 
GRS 1915+105 suggests a precessing warped ac- 
cretion disk, analogous to what is seen in other 
X-ray binaries. However, in all previous sources 
it was the soft X-ray flux which showed the long- 
term modulation. In contrast, in the case of GRS 
1915+105, it is the spectral slope of the power law 
and the radio flux. The radio flux is connected 
with the hard X-ray component and probably orig- 
inates near the inner part of the disk. It is unlikely 
that the periodic behavior of Fr and T can be ex- 
plained by a geometric effect of obscuration of the 
inner disk by the outer accretion disk. 

It has long been known (Katz 1973), that a pre- 
cession of a warped disk can be caused by tidal 
forces of the donor star. Assuming the accretion 
disk consists of small concentric rings which are 
tilted away from the orbital plane. Because of 
their fast rotation the rings behave gyroscopically. 
Without interaction of the various individual ring 
elements a retrograde precession around an axis 
perpendicular to the orbital plane will be excited 
by the tidal torque. The precession rate depends 
on the radius of the ring element which will lead 
to a tilted disk. In a fluid disk, an internal torque 
between the ring elements acts against the tilting. 
The resulting torque can enforce a precession of 
the disk. But the internal torque arises only if the 
disk is inclined with respect to the orbital plain 
and is connected with a dissipation of energy and 
a backward aligned movement of the disk into the 
orbital plane. 

For tidally forced precession, a correlation of 
Piong/Porb and q is theoretically expected but 
not found (Wijers & Pringle 1999). Neither the 
existence of a long term period nor the ratio 
Piong/Porb are correlated with the mass ratios of 
the components of the binary systems, ruling out 
tidal forces as the generally dominant mechanism 
for the long term periodicities in these systems. 
Thus, another mechanism is required which pro- 
vides continuous warping of the accretion disk. 
Radiation driven warping is one possible mecha- 
nism (Petterson 1977, Pringle 1996, Maloney et al. 
1996, Ogilvie & Dubus 2001). The outer, tilted 
part of the accretion disk is irradiated by the cen- 
tral source. If the radiation is absorbed and re- 
emitted parallel to the local disk gradient a torque 
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due to gas pressure affects the disk. Therefore, a 
small warp grows exponentially if the luminosity is 
sufficient to depress the dissipative processes forc- 
ing the disk back into the orbital plane (Petterson 
1977). 

In principle, for every system a combination 
of viscosity and accretion efficiency can be found 
which makes the system stable, unstable or highly 
unstable against radiation driven warping. A 
deeper knowledge of the internal physics of ac- 
cretion disk is therefore required to determine 
the cause of the warping. Recently, Ogilvie & 
Dubus (2001) found that radiation driven warp- 
ing in LMXB becomes relevant when the orbital 
period is above ~lday, which makes it possibly 
important for GRS 1915+105. 

Whether or not a warped accretion disk is the 
relevant cause of the long term periodicity in GRS 
1915+105 remains uncertain, because the 590 days 
period is visible in the non-thermal part of the X- 
ray spectrum and in the radio flux, but not in the 
disk parameters. However, it is worth mention- 
ing two issues: (i) a possibly existing long term 
variability in the soft X-ray flux and the disk pa- 
rameters can easily be masked by the high short 
term variability in timing and spectral behavior 
of GRS 1915+105. (ii) The observed periodicity 
can originate due to variability of the mass accre- 
tion rate or the viscosity of the inner part of the 
disk. One can assume that the Roche lobe over- 
flow of the donor star deposits it's mass at dif- 
ferent disk radii, depending on the mass transfer 
rate, thus differing mass density. The radial flow 
in the inner part of the accretion disk is therefore 
not constant, leading to a periodic variation of the 
amount of soft seed photons being comptonized in 
the corona. Therefore, the electron temperature 
in the corona varies on the same time scale due 
to Compton cooling, as the power law slope. The 
amount of matter to be ejected in a jet changes 
with similar periodicity which explains the overall 
periodicity in the radio flux. 

In the same way, periodic instabilities in the 
secondary affect the amount of accreted matter 
and may therefore cause or at least influence 
the observed periodicity. Wind-driven limit cy- 
cles (Shields etal. 1986) as discussed for LMC 
X-3 (Wilms etal. 2001) can also exist in GRS 
1915+105. In order to drive a wind from the 
outer parts of the accretion disk the sound speed 



for the inverse Compton temperature, kTjc, has 
to exceed the escape velocity in the outer disk. 
With an orbital separation of a ~ 7.5 • 10 12 cm 
the outer disk radii is likely of the order of ~ 
5 • 10 12 cm. Thus, GRS 1915+105 fuMls the cri- 
terium given in Equ.(3) of Wilms etal. (2001) for 
kTjc > 1 keV and may therefore be driving a sub- 
stantial Compton-heated wind. 

A warped accretion disk would temporarily 
shadow the irradiated donor star. A search for 
the 590 days periodicity in the IR flux of the donor 
would also help in the understanding of the origin 
and mechanism of the long term behavior of GRS 
1915+105. 

4. Conclusion 

We present the discovery of a 590 days long 
term periodicity in the power law slope of the hard 
X-ray spectrum, the hard X-ray flux and the ra- 
dio flux (at several frequencies) of the microquasar 
GRS 1915+105. No such periodicity is seen in the 
soft X-ray component. We discussed the observed 
behavior in the context of a warped accretion disk 
and of wind-driven limit cycles. We have found, 
that a periodic mass density variation in the disk, 
produced by the warp and/or by periodic instabil- 
ities in the donor, can describe the observations, 
by producing a varying amount of soft seed pho- 
tons for the Comptonization process and for the 
relativistic mass ejecta, very well. A wind-driven 
limit cycles is also a possible origin of the 590 days 
periodicity. 
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